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Abstract. The Laser Interferometer Space Antenna (LISA) guarantees the detection 
of gravitational waves by monitoring a handful of known nearby galactic binary 
systems, the so-called "verification binaries". We consider the most updated 
information on the source parameters for the thirty more promising verification 
binaries. We investigate which of them are indeed guaranteed sources for LISA and 
estimate the accuracy of the additional information that can be extracted during the 
mission. Our analysis considers the two independent Michelson outputs that can be 
synthesised from the LISA constellation, and we model the LISA transfer function using 
the rigid adiabatic approximation. We carry out extensive Monte Carlo simulations 
to explore the dependency of our results on unknown or poorly constrained source 
parameters. We find that four sources - RXJ0806.3+1527, V407 Vul, ES Cet and 
AM CVn - are clearly detectable in one year of observation; RXJ0806. 3+1527 should 
actually be observable in less than a week. For these sources LISA will also provide 
information on yet unknown parameters with an error between rs 1% and ~ 10%. 
Four additional binary systems - HP Lib, 4U 1820-30, WZ Sge and KPD 1930+2752 
- might also be marginally detectable. 



1. Introduction 

The Laser Interferometer Space Antenna (LISA) is an all-sky gravitational wave (GW) 
observatory in the frequency window 1CU 4 — 1 Hz with launch date in the 2015+ 
time frame. LISA is expected to observe gravitational radiation from a variety of 
sources P3 |2] , such as stellar mass binary systems and black hole binaries over a large 
mass spectrum. LISA guarantees the detection of GW's by monitoring several stellar 
mass (primarily white dwarfs) binary systems that are well known from electromagnetic 
observations and whose radiation is estimated to be sufficiently strong to be detected. 
They have been titled "verification binaries" and comprise systems from several classes: 
(i) mass-transferring AM CVn binary systems - consisting of a low mass donor star and 
a higher mass white dwarf accretor - in which mass transfer drives the orbital evolution 
to longer periods; (ii) double white dwarf binary systems, whose orbital evolution is 
primarily driven by gravitational radiation reaction, but possibly also tidal effects; (iii) 
(ultra-)compact X-ray binaries, in which a neutron star accretes material from a low 
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mass orbital companion; and (iv) cataclysmic variables (CV), consisting of an accreting 
white dwarf and (typically) a low-mass main-sequence star. 

Verification binaries are crucial because they guarantee the direct detection of 
GW's and provide a new opportunity to study the complex physics of compact objects. 
Previous analyses j3] of the detectability of verification binaries relied on crude estimates 
of their source parameters that are becoming progressively out of date due to new 
observations. As the mission formulation phase progresses, and firmer plans for data 
analysis are put in place, a more detailed study is needed to identify those systems that 
can be indeed considered "verification binaries" (i.e. guaranteed sources) and to explore 
what kind of new information LISA will be able to provide. 

In this paper we consider the most updated information on the thirty most 
promising nearby galactic binaries; we carry out a systematic study of the signal- 
to-noise ratio (SNR) at which they can be observed and explore the accuracy with 
which LISA can measure the unknown source parameters, keeping into account prior 
information. Our analysis shows that four systems (RXJ0806. 3+1527, V407 Vul, ES Cet 
and AM CVn) can indeed be considered guaranteed sources and, as such, verification 
binaries. Four others (HP Lib, 4U 1820-30, WZ Sge and KPD 1930+2752) might also 
be marginally detectable; this will depend on the actual polarisation of gravitational 
waves, the duration of the mission, the exact noise level and the performance of data 
analysis algorithms. We also show that for RXJ0806.3+1527, V407 Vul, ES Cet and AM 
CVn, LISA will provide new information on unknown parameters (such as the orbital 
inclination angle and a combination of the distance and chirp mass) with an error ~ 1% 

- 10%; this will be important to shed new light on the physics at work in such compact 
objects. The remaining binaries appear to be beyond the reach of LISA. 

2. Sources and signal model 

A large number of nearby galactic compact binary systems have been discovered and 
monitored by recent surveys; 30 of them have been identified as promising (or certain) 
candidates for GW detection with LISA: we will refer to them as potential verification 
binaries. In this section we introduce the sources and the relevant information that are 
currently available. We also review the model of the GW signal that is used in the next 
section to study the expected SNR characterising LISA observations and the estimates 
of the statistical errors associated to the measurement of unknown parameters. 

In Table we summarise the 30 most promising sources for LISA and their 
parameters: more information, including references can be found in []\ . The location of 
the systems in the sky - identified by the unit vector N and ecliptic coordinates (9^, (fx) 

- and the orbital period (and therefore the frequency / of GWsif) are extremely well 
determined; they can be considered (at least for the analysis presented in this paper) 
exactly known. On the other hand, some parameters that are essential to evaluate the 

I GW's are emitted at all the multiples of the orbital frequency: here we just consider the leading 
order Newtonian mass quadrupole radiation, that provides by far the dominant contribution. 
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Coordinates 


/ [mHz] 


D [pc] 


m 2 [Mq] 


mi [M e ] 


AM CVn systems 


RX J0806. 3+1527 


120.443557 


-4.704035 


6.2202766 


300-1000 


0.13 


0.2-0.5^1 


V407 Vul 


294.994600 


+46.783096 


3.51250 


300-1000 


0.068 


0.70 


ES Cet 


25.310039 


-20.333218 


3.22 


350-1000 


0.062 


0.70 


AM CVn 


171.084432 


+37.441925 


1.94414 


606 


0.14 


0.85 


HP Lib 


235.786580 


+4.959494 


1.813 


197 


0.032 


0.57 


CR Boo 


202.971814 


+17.896456 


1.360 


337 


0.023 


0.55 


KL Dra 


334.830867 


+78.322322 


1.333 


100* 


0.022 


0.27 


V803 Cen 


216.866251 


-30.317527 


1.241 


347 


0.021 


1.31 


SDSS J0926+3624 


132.286781 


+20.234177 


1.177 


100* 


0.02* 


0.6* 


CP Eri 


42.830979 


-26.426918 


1.176 


100* 


0.019 


0.63 


2003aw 


140.733734 


-21.234213 


0.9862 


100* 


0.015* 


0.42* 


SDSS J1240-0159 


190.193388 


+2.225887 


0.8921 


350-440 


0.015* 


0.38* 


GP Com 


188.418156 


+23.000197 


0.7158 


75 


0.010 


0.45 


CE 315 


206.451728 


-14.462612 


0.5120 


77 


0.006 


0.48 



(Ultra-)compact X-ray binaries 



4U 1820-30 
4U 1543-624 
4U 1850-087 
4U 1626-67 
CC Com 


275.841980 -7.0273940 
251.159958 -41.352944 
283.530386 +14.112497 
259.038082 -44.908517 
174.038625 +21.7757800 


2.920 
1.832 
1.618 
0.4GS1 

0.105 


81000 
5000* 
8200 

80000 
900 


<O.10 
0.04 
0.03 
0.02-0.08GS1 
0.36E1 


1.40 

1.4 

1.4 
I4EEI 
0.62GH 


Double white dwarfs 


WD 0957-666 
KPD0422+4521 
KPD1930+2752 
WD 1101+364 
WD 1704+481 
WD 2331+290 


209.226640 -67.301935 
76.322155 +22.854567E1 
302.454582 +48.9048890 
152.950020 +27.689153 
242.213221 +70.2247280 
7.631717 +29.202332 El 


0.379520267 
0.25669oE] 
0.2434262GH 

0.15996 
0.1598789E3 

0.1 EH 


135 
100* 
100* 

97 
100* 
100* 


0.32 
0.51lGSl 
0.5EI 

0.31EI 
0.39G! 

>0.32 GH 


0.37 
0.526GS1 
0.97G1 

o.seCni 

0.56G1 
0.39G5I 


Cataclysmic variables 


EI Psc 

SDSS J1507+5230 
GW Lib 
WZ Sge 

SDSS J0903+3300 


356.370891 +8.925019 
192.423733 +64.768925 
234.794265 -6.417916 
309.727920 +36.921088 
128.734149 +15.542558 


0.5195 
0.4958 
0.4341 
0.4065 
0.3918 


210 
100* 
100* 
43 
100* 


0.13 
0.13* 
0.13* 
<0.11 
0.13* 


0.7 
0.7* 
0.7* 
>0.7 
0.7* 



Table 1. The 30 close stellar mass binary systems that are currently regarded 
as potential "verification binaries" for LISA. The table summarises the relevant 
parameters for LISA observations: ecliptic coordinates (J2000.0), gravitational wave 
frequency /, distance D to the source and masses mi and vn-i of the two stars. The table 
is based on the observational summary reported in [2] and provides explicit references 
where additional information have been included. For those systems whose distance 
is not known we set D — 100 pc, a lower limit to D due to the lack of detection of 
proper motion. For 4U 1543-624 we adopt D = 5000 pc 0] due to its astrophysical 
nature. These values are marked by a superscript star. In six cases - SDSS J0926, 
2003aw, SDSS J1240, SDDS J1507+5230, GW Lib and SDDS J0903+3300 - the mass 
determination is highly uncertain or not available. For them (again marked by a 
superscript star) we adopt canonical values for this class of systems and keep into 
account estimates of the mass ratio (see the text for more details). In the table we 
classify RX J0806.3+1527 and V407 Vul as AM CVn systems, but we note that their 
classification is still debated, see 0] and references therein. 
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expected SNR are either unknown or poorly constrained. For none of the sources we have 
information about the orientation of the binary orbital angular momentum L (described 
by the polar coordinates 9l and <Pl), which determines the signal polarisation, and in 
turn the actual GW amplitude at the detector through the antenna beam patterns. 
For several systems we also have only poor (if any) information about the distance 
to the source D and/or the masses of the stars mi and m 2 , which affect the GW 
amplitude. In this paper, we adopt a distance D = 100 pc for KL Dra, SDSS J0926, CP 
Eri, 2003aw, KPD 0422+4521, KPD 1930+2752, WD 1704+481, WD 2331+290,SDDS 
J1507+5230, GW Lib and SDDS J0903+3300; note that this value represents a lower 
limit to D, because the proper motion of those objects is not observed. For 4U 1543-624 
we adopt a distance D = 5 kpc j3], based on evolutionary arguments. For SDDS J0926 
no information on the masses are available: in our analysis we assume typical mass 
values for AM CVn binaries [IZj, i.e. m 2 = 0.02 M® and mi = 0.6 M . For 2003aw and 
SDDS J1240, some measurements of the mass ratio have been obtained (r^/mi ~ 0.036 
and ~ 0.039 for 2003aw and SDSS J1240, respectively): by keeping into account these 
constraints and using a standard AM CVn model jTTJ we adopt the mass values as 
reported in Tabled For the three CV systems SDDS J1507+5230, GW Lib and SDDS 
J0903+3300, whose masses have not been measured, we set the values of mi and m 2 to 
those of EI Psc, the only CV with a reliable mass determination. 

In our analysis we consider the output of the two uncorrelated Michelson channels 
as introduced in JH] and include the effect of the LISA transfer function using the rigid 
adiabatic approximation ^H]- The gravitational wave signal h(t; A) at each of the two 
readouts can therefore be written as 



where 



h{t;\) = '£B n (t)cosXn(t), (1) 



n=l 



B n (t) = [(F^\t)A + (t)f + (F^(t)A x (t)r] 1/2 , (2) 

Xn(t)=(f> gw (t) + Mt) + Mt)- (3) 

The Doppler phase shift 0d(^) is caused by the motion of the interferometer around the 
sun and the terms 

F^{t)A x {t) 



4> n (i) = arctan 



(4) 



. Fi + \t)A + (t) 

are the polarisation phases induced by the change of LISA's orientation. i^ + ' x - ) (t) are 
the antenna beam patters. Note that 0„(t) and F^ + ' x \t) are different for the two 
Michelson channels. The amplitudes of the two independent polarisations in Eq. (j2J) are 
given by 

A+it) = 1 J — [i+(L-N)\ [rrf(t)f 3 , (5) 
A/f 5 / 3 

A*(t)= -4^- (L-N) [nf(t)f 3 , (6) 
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where the chirp mass is defined as M. = (mi + m 2 ) 2//5 [(m 1 m 2 )/(m 1 + m 2 )] 3 / 5 . The 
amplitudes depend on the inclination angle i given by 

cos i = L ■ N = cos 9l cos On + sin Ol sin On cos(y?L — yjjy) • (7) 

In this paper we model the gravitational wave signal from the sources in Tabled as 
exactly monochromatic with respect to an observer at rest in the solar system barycentre. 
In reality the frequency does change due to radiation reaction, mass transfer and tidal 
effects. However, we estimate that for the vast majority of the sources considered here 
the frequency drift will not be observable during the LISA mission. Exceptions are likely 
to be RXJ0806. 3+1527 and V407 Vul, where measurements of the period derivatives 
have already been reported jlj; in fact, LISA could provide new direct measurements 
of their orbital evolution and details on the physics at work. Neglecting this intrinsic 
frequency drift (as we do in the present analysis) does not introduce any sizable effect 
on the SNR that is computed in the next section, and it affects the estimates of 
the statistical errors of the measurements only in a negligible way compared to the 
present uncertainties on the source parameters, see [2*T] . In summary the phase of the 
gravitational wave signal in Eq. (J3J) is given by 

gw (t) = 2nf t + $ , (8) 

where fo and $o are the frequency and phase, respectively, at the beginning of the 
observations. It is also useful to introduce the constant amplitude 

on M 5/3 

A = 2(Kf ) 2/3J —. (9) 

The signal h(t; A) that we consider is hence described by the parameter vector 

A = {ha.A,^Q,fo,cosO N ,ip N ,cosO L ,ip L }. (10) 
We introduce the usual inner product between two signals h and g as 



r-+oo h*( f)a( f) 2 
(%) = 2/ { : !Z df ~ — / h(t)g(t)dt, (11) 



oo 



Sn(f) So 

where the latter equality follows from Parseval's theorem; S n (f) is the one-sided noise 
power spectral density that we consider constant around / ; we introduce the notation 
So = S n (fo)- In our analysis we model the noise spectral density according to the 
latest estimate of the instrumental and "confusion noise" [22]- in the noise produced 
by astrophysical foregrounds we include unresolved galactic and extra-galactic white 
dwarf binaries, but not the contribution from captures of solar-mass compact objects by 
massive black holes. If such foreground radiation is indeed present at the level described 
in [22], the SNR at which LISA can detect the verification binaries with /o in the mHz 
region would be reduced by a factor ^ 1.5 with respect to what we report in the next 
section. The optimal signal-to-noise ratio p at which h(t; A) can be detected is given by: 

p 2 = ]>> {a W a) ) , (12) 



a=l 
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Figure 1. The signal-to-noise ratio at which LISA can observe the known binaries in 
one year of observation. We show the combined optimal signal-to-noise ratio p from the 
two independent LISA Michelson outputs. The diamonds correspond to the median 
of the distribution of p over a population of 10 4 sources, whose parameters have been 
chosen according to Tabled (see text for more details). Open diamonds are used for 
those systems for which we have adopted "fiducial" values for the distance and/or 
chirp mass (the parameters labelled by a superscript star in Table The error bar 
is drawn to cover the range between the minimum and maximum value of p obtained 
over the 10 4 Monte Carlo trials. The right panel refers to the potential verification 
binaries that have been classified as AM CVn, the left panel to the remaining systems. 



where a = 1, 2 labels the Michelson outputs. The statistical error AA associated to the 
measurement of A can be computed using the variance-covariance matrix. In the limit 
/9> 1, A A follows the Gaussian probability distribution 



1/2 



where the Fisher information matrix is given by 



a \ 



dX k 



(13) 



(14) 



The (lower limit to the) expected mean square error on each parameter can be computed 

as 



<(AA^) 2 >= (r) 



J.I 



(15) 



Note that the errors depend on the actual value of A and scale with the SNR according 
to ((AX j ) 2 ) oc 1/p 2 . 
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3. Results 

In this section we report the results of the exploration of the SNR and the expected 
mean square errors on the parameter measurements for the sources in Table [0 Due to 
the uncertainty on some of the parameters that determine the actual signal received at 
the detector - L for all the sources and, in a number of cases, D and M. - we performed 
Monte Carlo simulations to explore the range of the values that can be obtained during 
the actual LISA mission. For each binary we consider 10 4 realisations of the same 
source with fixed position N and GW frequency f (according to the values reported 
in Table Q). For each realisation, L is chosen randomly and drawn from a uniform 
distribution over —1 < cosOl < +1 and < ipi < 2n. In addition, for those systems 
for which D is known within a given range (RX J0806. 3+1527, V407 Vul, ES Cet and 
SDSS J1240) the distance is randomly selected from a uniform distribution within the 
relevant interval. The same is true for mi for RX J0806. 3+1527. For three sources 
(4U 1820-30, WD 2331+290 and WZ Sge), only limits to m x and/or m 2 are available; 
in all three cases we fix in the Monte Carlo simulations the value of rrii and m 2 to the 
actual limits reported in Tabled Finally, we use for the systems with partially unknown 
parameters the assumed fiducial values for D, mi and m2, as marked by a superscript 
star in Table H We assume the standard observation time T = 1 yr, and for each of 
the realisations we compute the optimal SNR, Eq. (|T2"|) . and the errors on parameter 
extraction, Eq. (JT3J). For the latter, we always assume f and N to be perfectly known; 
we therefore compute a four dimensional Fisher information matrix for A, <£>o, @l and 
ifL. We also evaluate the error on the inclination angle l, cf Eq. (J7j). In summary for 
each verification binary we have 10 4 values of SNR and error on each of the parameters: 
we then compute the median over those values to obtain a "typical" result. For the 
SNR we also record the minimum and maximum value and use this as a measure of the 
spread in the results. For AA this is not possible because for some choices of L some 
of the parameters become degenerate (and we obtain unphysically large values for AA). 
Instead we compute the sample standard deviation of AA- 7 (for each j = 1,2,3,4) and use 
it as error bar on the typical values. The results are summarised in Figures 1, 2 and 3. 

We consider first the SNR in order to determine which of the potential verification 
binaries will be clearly observable. The results are shown in Figure ^ The exact SNR 
threshold to claim detection depends on a number of details that go beyond the scope 
of this paper. However, as the search parameters are almost exactly known (source 
position and gravitational wave frequency), we conservatively set it to p = 5. If we 
consider the median value of the SNR, p, as the figure of merit to be compared to the 
threshold, then only four AM CVn systems are detectable: RXJ0806. 3+1527 (p = 62), 
V407 Vul (p = 30), ES Cet (p = 19) and AM CVn (p = 8). If the orientation of L is 
favourable, LISA could reach a value of p as high as 227, 79, 62 and 13, respectively. 
It is important to notice that RXJ0806. 3+1527 will likely be detected within a week 
and V407 Vul and ES Cet within a month of science operations. We consider these 
results robust, because these four systems emit in the frequency band 3 mHz - 7 mHz, 
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Figure 2. The expected minimum mean square error that characterises the 
determination of the GW amplitude of the known binaries using one year of LISA 
data. We plot the relative error AA/A - the amplitude A is defined according to 
Eq. © - as a function of the GW emission frequency. The diamonds correspond 
to the median of the distribution of AA/A over a population of 10 4 sources, whose 
parameters have been chosen according to Table H (see text for more details) . Open 
diamonds are used for those systems for which we have adopted "fiducial" values for the 
distance and/or chirp mass (the parameters labelled by a superscript star in Table 
The error bars correspond to the sample standard deviation as measured from the 
10 4 Monte Carlo trials. The right panel refers to the potential verification binaries that 
have been classified as AM CVn, the left panel to the remaining systems. 



in which sources are sufficiently sparse and a possible additional contribution to the 
noise from confusion as generated by extreme-mass ratio inspirals (if indeed present) 
degrades the SNR by less than a factor of 2 |22j- Furthermore, the use of the Time 
Delay Interferometry combinations tuned to a given source should increase the SNR by 
a few tens of percent, see f2Hj and references therein. We also note that one AM CVn 
(HP Lib), one LMXB (4U 1820-30), the Cataclysmic Variable WZ Sge and the double 
degenerate white dwarf KPD 1930+2752 would yield a maximum SNR 5, but with 
the relevant median values below 3. Given the GW emission frequency of these systems, 
a marginal detection depends on details - such as degree of signals overlap and actual 
performance of data analysis algorithms - that will be known only at the time of the 
real analysis. 

The remaining sources that have been so far included into the potential verification 
binaries will likely not be observable by LISA; however, several parameters are 
sufficiently uncertain that a conclusive statement in this respect is premature. We 
investigated also the likelihood of observing the double radio pulsar PSR J0737-3039: 
unfortunately, even for a 10 years mission the SNR is only w 0.3, which suggests that 
this unique object will be outside LISA's ken. 
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Figure 3. The expected minimum mean square error that characterises the 
measurement of the inclination angle of the known binaries using one year of LISA 
data. We plot the error A cos t, where cos l is given by Eq. , as a function of the 
GW emission frequency. The symbols are the same as those in Figure 

LISA will provide a direct measure of the amplitude A, see Eq. ©, thus of the 
ratio Ai 5 / 3 /D, and of the elusive (for electromagnetic observations) inclination angle l, 
offering a new opportunity to study the physics of compact objects and white dwarfs 
in particular. Furthermore, if either the chirp mass or the distance are independently 
known, one could extract information on the other parameter from A. The possibility of 
measuring M. is particularly important for AM CVn systems because it could shed new 
light on the physical effects responsible for the binary orbital evolution [23]. in Figures^] 
and 01 we show the range of the expected mean square errors A A/ A and A cost. For 
the clearly detectable verification binaries we have 0.02 ^ AA/A ^ 0.3. Assuming that 
independent measurements of the distance are of comparable quality by the time LISA 
is in operation, one would be able to extract information on the chirp mass with a ~ 10% 
error. The error on cos l is also sufficiently small, suggesting that one could measure 
cos i with an error in the range m 0.01 — 0.1. 

4. Conclusions 

We have shown that there are indeed four (possibly eight) verification binaries for LISA. 
They will be surely detected within the first year of the mission, with one of them, 
RXJ0806. 3+1527, observable by the end of the first week of science operation. High 
quality observations will be possible for RXJ0806. 3+1527, V407 Vul, ES Cet and, to 
a lesser extent, AM CVn. Our results are clearly affected by uncertainties regarding 
the actual level of noise during the mission and in particular the effect of astrophysical 
foregrounds. Furthermore, in our analysis we have not included the full power of the 



The LISA verification binaries 



10 



technique of Time Delay Interferometry which will be crucial for the real analysis of the 
data and will improve the results that we have quoted here by a few tens of percent. 
Nevertheless, we consider the main conclusions of our paper to be robust. It is now 
important to develop an end-to-end analysis strategy that can deliver the expected 
results on real data. This work is currently in progress and will be reported elsewhere. 
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